November, 2006 



Vacuum Birefringence in a 
Rotating Magnetic Field 

Stephen L. Adler 
Institute for Advanced Study 
Princeton, NJ 08540 



Send correspondence to: 

Stephen L. Adler 

Institute for Advanced Study 

Einstein Drive, Princeton, NJ 08540 

Phone 609-734-8051; FAX 609-924-8399; email adler@ias.edu 



1 



ABSTRACT 

We calculate the vacuum polarization-induced ellipticity acquired by 
a linearly polarized laser beam of angular frequency lD on traversing a region 
containing a transverse magnetic field rotating with a small angular velocity 
Q, around the beam axis. The transmitted beam contains the fundamental 
frequency Co and weak sidebands of frequency d)±2Q, but no other sidebands. 
To first order in small quantities, the ellipticity acquired by the transmitted 
beam is independent of Q, and is the same as would be calculated in the 
approximation of regarding the magnetic field as fixed at its instantaneous 
angular orientation, using the standard vacuum birefringence formulas for a 
static magnetic field. Also to first order, there is no rotation of the polar- 
ization plane of the transmitted beam. Analogous statements hold when the 
magnetic field strength is slowly varying in time. 
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1. Introduction 

There has recently been renewed interest in trying to measure the birefringence of the 
vacuum in a strong magnetic field predicted by the Heisenberg-Euler effective Lagrangian 
[1,2,3,4] for quantum electrodynamics (QED). Such a measurement was a primary motivation 
for the PVLAS experiment [5], which however has reported a signal about 10^ times larger 
than expected from QED. In a recent article [6], a claim was made that just the kinematics 
of rotation of the magnetic field used in the PVLAS experiment, through the generation of 
multiple sidebands to the transmitted laser beam, could explain the PVLAS result. However, 
examination of ref [6] and a related paper [7] suggests several problems [8]. First, the authors 
have overestimated the strength of the vacuum birefringence coefficient by a factor of over 
1000. Second, since on bases rotating with the magnetic field, the vacuum as modified 
by the magnetic field still acts on an incident electromagnetic wave as a time-independent 
linear medium, the presence of an infinite sequence of sidebands (as would be obtained from 
a nonlinear medium) is not expected; the only sidebands should be those associated with 
transformations between fixed laboratory bases and the rotating bases. 

Nonetheless, the authors of ref [6] have called attention to an interesting, and as it 
turns out quite nontrivial, problem in wave propagation. This paper is devoted to formulating 
and solving this problem to leading order, using standard scattering theory and perturbation 
theory methods. Our conclusion is that to leading order in small quantities, the method used 
by the PVLAS group to calculate the birefringence-induced ellipticity yields the correct 
answer, which is independent of the rotation angular velocity Q. 

To set the stage, we start from the lowest order formulas describing electromagnetic 
wave propagation in an external magnetic field, as deduced from the Heisenberg-Euler ef- 

3 



fective Lagrangian. With repeated indices summed, the d and h fields are obtained from the 
e and b fields by 

ds —^st^t ) 

(la) 

The polarization tensors e^t and //^^ arising from the fourth order box diagram are given by 

(16) 

/.;/=5,,(i-20-4e4(t)A(t) , 

with Bg{t) a unit vector along the magnetic field, which we take to have fixed magnitude B. 
With B measured in unrationalized units, and setting h — c — 1, the parameter ^ governing 
the strength of the QED vacuum polarization effect is given by 

where a ~ 1/137.04 is the fine structure constant and m is the electron mass. Finally, in 
terms of d, h, e, and b the Maxwell equations take the form 

V-d^V-b^O , 

^ - db ^ r 9d ^^^^ 
Vxe = -^, VxA = - . 

Equations (la-d) are the fundamental physical equations governing the wave propagation 
problem under study. When the magnetic field orientation B{t) is time-independent, anal- 
ysis of this problem shows that the vacuum becomes birefringent. Specializing to the case 
of a wave-vector k perpendicular to B, the polarization state with e perpendicular to B 
propagates with index of refraction k/cu — 1-1-2^, while the polarization state with e parallel 
to B propagates with index of refraction k/u = 1 + |{ (with k = \k\). We will denote the 
difference between the larger and smaller refractive indices by An = 



To finish specifying the problem, we shall assume that the linearly polarized laser 
beam is propagating upwards along the z axis, with field strengths (expressed in terms of 
(x, y, z) components) given by 

e =d = e^'^("-*)(cos^,sin^,0) , 

(2a) 

h^h^ e^'^(^-*) (- sin cos 0) . 
We assume that the the uniform magnetic field is confined to the region < 2; < L, and is 

oriented perpendicular to the z axis, and rotates uniformly with angular velocity according 

to 

B (t) = (cos , sin nt,0) , (26) 

from which we learn that (with the prime denoting time differentiation) 

B'{t) ^n{-sinnt,cosnt,0) ^nz X B . (2c) 

The problem is then to calculate the upward moving transmitted wave emerging at ^ = L. 

Before proceeding with details, we outhne our basic strategy. The first observation 
is that since Eq. (2c) relates the time derivative of B to z x B with a time-independent 
coefficient, the wave propagation problem in the region < z < L will involve differential 
equations with time-independent coefficients when referred to the orthonormal rotating bases 
B and z x B. These differential equations can then be solved by a standard traveling wave 
Ansatz; this calculation of the propagation modes in the rotating magnetic field "medium" is 
carried out in Sec. 2. One then has to transform the incident wave to the rotating bases, and 
do a matching of incident, transmitted, and reflected waves at the medium boundaries z — 
and z = L. This calculation is carried out in Sec. 3. The calculation is greatly facilitated 
by working to leading order in the small quantities ^, Q/o), and QL, that is, we regard both 
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^ and Q as small. Since each reflected wave is reduced in strength by O(^), this allows us 
to neglect multiple reflections, such as a wave reflected back from z = L, and then reflected 
forward from z — 0, and finally emerging ai z — L, which makes a contribution of order 
to the transmitted wave. When multiple reflections are neglected, one can solve the wave 
propagation problem (to leading order) by considering independent matching problems at 
z — and z — L. Moreover, it then suffices to compute the phases of the various waves to 
first order in small quantities, but the corresponding transmitted wave amplitudes only to 
zeroth order. The reason this approximation is adequate is that the transmission coefficients 
at the two boundaries are independent of L when multiple reflections are neglected, and so 
first order terms in the transmission coefficients can contribute only an L-independent term 
to the polarization parameters of the emerging wave ai z = L. However, when L=0 no 
birefringent medium is traversed, and so the ellipticity of the emerging wave is zero; hence 
the first order corrections to the wave coefficients must make a vanishing contribution to 
the ellipticity, and similarly, to the rotation of the polarization axis. The calculation of the 
ellipticity and polarization axis, from our formulas for the emerging waves, is given in Sec. 4. 
An alternative calculation method, based on a perturbative expansion in ^, is given in Sec. 
5, which gives results identical to the wave- matching calculation of Sec. 3 for the rotating 
field case, and also applies to the case when the magnetic field strength is time dependent. 
A brief discussion of our results and their experimental implications is given in Sec. 6. 

2. Electromagnetic wave eigenmodes in a rotating magnetic field 

We begin by calculating the electromagnetic wave eigenmodes in the region contain- 

— * 

ing the rotating magnetic field. We write the e and b fields in terms of components along 
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the rotating bases B and z x B, 

e=[Ei{t)B + E2{t)z X B^' , 

(3a) 

b^[Bi{t)B + B2{t)zx B]e'''' . 
Substituting these into Eq. (la), and using Eq. (lb), we get the corresponding expressions 

for the d and h fields, 

d=[{l + 5i)Ei{t)B + (1 - 2C)E2{t)z X B]e"'' , 

(36) 

h =[(1 - 60Bi{t)B + (1 - 2i)B2{t)z X BY^' . 
Substituting Eqs. (3a,b) into the MaxweU equations of Eq. (Id), and using Eq. (2c), we get 

a set of coupled equations for the coefficient functions Ei^2{t) and Bi^2{t), 

ikEiit) = - VtBi{t) - B'^{t) , 

-ikE2it)^-B[{t)+nB2it) . 

(4a) 

ik{l - 60B,{t) =(1 + 50^E,{t) + (1 - 20E'^{t) , 

-zk{l - 20B2{t) =(1 + 50E[{t) - (1 - 20nE2{t) . 
Since this set of equations has time independent coefficients, it can be solved by making an 

exponential Ansatz, 

^1,2 (t) = ^1,26-^"* , B,,2{t) = B,,2e-'^' . (46) 

Substituting Eq. (4b) into Eq. (4a), and simplifying by neglecting second order small terms 
proportional to and ^Jl, we get the coupled equations for the constant amplitudes £'1^2 
and Bi^2, 

kEi ^ujB2 + niBi , 
k iE2 = — uriBi — VLB2 , 

(4c) 

kiBi^- u{l + 40 iE2 + ^Ei , 

kB2^uj{l + 7^)Ei-niE2 . 
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We note at this point that we have ignored effects arising from the electric field induced 
by the rotating magnetic field; these are of order ^fl at least, and so are second order 
in the small parameters of the problem. Similarly, other relativistic effects associated with 
electrodynamics in rotating frames, on which there is a substantial literature, are not relevant 
for our analysis. 

Since Eq. (4c) is a homogeneous set of linear equations, a solution is possible only 
when the determinant vanishes, which yields a quartic equation giving the dispersion relation 
of the wave, 

[k'^ -u}^{l + 70][k^ -Lj\l + 40]-2fl^{uj^ + k^) + n^ ^0 . (5a) 

This can be solved to give two solutions k± which describe upward propagating waves, and 
two solutions —k± which describe downward propagating waves, with k± given by 

k± —UJ ± Q + -0J(T± , 

, (5&) 
a± =(7±(u;) = y e ± W^/^f + (3^/2)') ^ - , 

so that the corresponding refractive indices are given by 

k±/uj — n± — 1± h -(7± . (5c) 

Note that (t± is first order in small quantities, and is uniformly of order ^ or smaller irrespec- 
tive of the relative sizes of ^ and VL/uj. Also, we note that the terms of order ^VL that we have 
dropped make a contribution inside the square root of order ^(VL/uY (because the dispersion 
relation must be an even function of Q, as we have verified by explicit computation), and 
so is a higher order effect relative to the quadratic terms that we retained inside the square 
root. 



Solving for the ratios of the field coefficients to Ei, we get for the case of upward 
propagating waves with oscillatory behavior e'("±^~*)'^ 

— =^1 + a± , (6a) 

with the quantities a±, (3±, j±, which all vanish when ^ vanishes, given by 

- - a± , 

where the first lines of the expressions for a± and 7± are accurate to first order in small 
quantities, and the second lines are accurate to leading zeroth order. For the case of down- 
ward propagating waves with oscillatory behavior e'(~"±^~*)'^, one simply reverses the sign 
of B2/E1 and iBi/Ei, while keeping iE2/Ei the same as in Eqs. (6a,b). This completes the 
analysis of the electromagnetic wave eigenmodes in the region < 2; < L containing the 
rotating magnetic field B. 



3. Wave matching to get the transmission coefficient 

Now that we know the form of propagating waves in the magnetic field region, we 
turn to the problem of determining the emerging waves at 2; > L produced by the laser beam 
incident from below z = 0. As noted in Sec. 1, when doubly reflected waves are ignored, 
we can solve this problem by treating sequentially the wave matching at 2; = 0, followed by 



the wave matching ai z — L. The first step is to rewrite the incident wave at z = on the 
rotating basis, giving 

e1,=o =e-^^*[cos(^ -nt)B + sm{e -Vtt)zxB] , 

(7) 

&|^=o ^e-'^^- sm{d - nt)B + cos{9 - nt)z x B] , 

which shows that all the components of the incident wave on the rotating bases have the 
time dependence e~*('^+^)* or e"'^'^"^)*. Thus to do the match, we need a reflected wave at 
2; < with four coefficients (two for the refiected wave with frequency cu + ft, and two for 
the reflected wave with frequency lj — fl), and a transmitted wave with four coefficients (one 
each for the ± eigenmodes with frequency cu — cu + D,, and one each for the ± eigenmodes 
with frequency lu = cu — fl). Continuity of the e and h fields across the junction at 2; = 
gives eight matching conditions (there are two field components parallel to the boundary 
for each of the two fields to be matched, and two frequency components for each). Solving 
the wave matching equations shows, as expected, that the refiected wave amplitudes are of 
order ^, so to get the transmitted wave amplitudes to zeroth order it suffices to ignore the 
reflected waves. Similarly, a.t z = L there are also eight unknown coefficients, four for the 
transmitted wave (two for frequency a) + Q and two for frequency a) — Jl) , and four for the 
waves of the two frequencies reflected back into the magnetic fleld region z < L. Again there 
are eight matching conditions, and solving for the reflected waves shows that they are again 
of order ^ in magnitude, and so again the leading order coefficients of the transmitted wave 
can be obtained by ignoring the refiected waves. 

Thus, the transmitted wave on the rotating basis, to leading order in the wave am- 
plitudes, is simply the result of compounding transmission at z = 0, followed by propagation 
from z — to z — L, and then followed by transmission at z — L. The final step is to 
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convert back from the rotating bases to fixed laboratory bases. This yields three distinct 
emerging wave components, with frequencies uj, uj + 2VL and uj — 2Vt. We dispense with all 
of the intermediate algebra, and simply present the final result for the transmitted wave at 
z = L, which takes the form 

&l.=L=e^''^^-*)(-y,x,o) , 

with the components X, Y given in terms of quantities A,B,C,D by 

1 



X ^-{C + D + Ae^'^^ + Be-^'^*) , 

Y ^-(-C + D - Ae'''^' + Be-^'^') . 
2 



Introducing the definitions 



(8a) 



(86) 



(8c) 



(9a) 



the auxiliary quantities A,B,C,D are given by the following formulas, 

^4 -(^c^'^+L ^i(-2n+^-)L^^^-ie (1 + Q-)q+ 

1 + a_ — a+ ' 

g ^(^i(2n+4>+)L (.i't'-LyW (1 ~ Q+)Q- 

1 + a_ — a+ ' 

(J ^^ie+i(2n+ct>+)L a+a- _^ ^ie+i4>-L ~ + Q-) 

1 + a_ — a+ 1 + a_ — a+ ' 

1 + a_ — a+ 1 + a_ — a+ ' 

These formulas are accurate to first order in small quantities in the phases, and to zeroth 
order in small quantities in the real amplitudes multiplying the phases. As we stressed in 
Sec. 1, first order corrections to the real amplitudes are independent of L, and so can make no 
contribution to physical attributes, such as the ellipticity, of the emerging wave. Expanding 
the exponentials through first order in small quantities, we get the following expressions for 
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A, B, C, D (with (7± = a±{{D)) , 
A —icuL 



2n 1, 



B —iujL 
C 
D 



2Q 1. 



(1 + a_)a+ 
1 + a_ — a+ 
(1 — a+)a_ 



, 2n 1 

* (D 2 +y 1 + a 



1 + a_ — a+ 
a+a_ 



. 1 (l-a+)(l + a_) 

+ lujL-a-—^ — — 

2 1 + a_ — a+ 



(96) 



^ . 1 (l-a+)(l + a_) . / 2(] 1 

2 l + a_ — a+ \ ^ ^ /l + a_ — a+ 



These equations can be greatly simphfied by using the following three algebraic identities 

obeyed by the quantities a±, 

a+ — a_ + 2a_|_a_ =0 , 



^1 — a+)(l + a_ 
1 + a_ — a+ 



-1 - 



1 + a_ 



(10a) 



1 + a_ - a+ =1 + ^(c^+ - 0--) 



We then get the remarkably simple formulas 



11 3 

X^cose{l + — icD-^O + -iu}L^cos{29.t - 9) 

11 3 
r =sin6'(l + —iOL^) + -iuL^sm{2nt - 9) 



(106) 



which together with Eq. (8a) are our final result for the wave transmission problem. They 
will be used in the next section to calculate the ellipticity and major axis orientation of the 
emerging wave. 

The vacuum with magnetic field is a specific instance of a generic weakly birefringent 
medium which has orthogonal polarization eigenmodes with refractive indices n\\, n±. Hence, 
with the replacements — > ny — 1, 2^ ^ n± — 1, — > + n± — 2, and |^ — >• ny — n±, 
the preceding discussion generalizes to the case in which such a generic medium is slowly 
rotated. 
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4. EUipticity and polarization cixis of the emerging wave 

We proceed next to calculate the polarization parameters of the emerging wave, 
following the exposition of Born and Wolf [9] . Let us write 

X =cose + iA^x - cos 66'^'^/''°'^ = aie'^' , 

(11a) 

Y = sin 9iA^ ~ sin ^e^^^/^'"^ = ose^*^ , 
with A^ y first order small quantities that can be read off from Eq. (10b). Then the auxiliary 
angle a of Born and Wolf is given by 

tana = — = tan^ , (11^) 

which implies that a = ^ up to second order corrections, and so to first order there is no 
rotation of the polarization axis of the beam. The second auxiliary angle x oi Born and 
Wolf is given by sin 2% = sin 2a sin S; since 5 — 82 — Si — Ay / sin ^ — A^ / cos 9 is first order 
small, to leading order (and allowing the ellipticity to carry a sign) we have 



^ =Ellipticity = tanx — x — -5sin 2q; ~ -5 sin 2^ 

= cos 61 Af. - sin 61 A^ . 
From the values of A^y obtained from Eq. (10b), we then find 



(11c) 



^ ^ ~^CuLsm2{e -Qt) , (lid) 

which using Cu — 27r/A and ^ = 2 An/3 is equivalent to 

ttL 

* = — -Ansm2{0 -fit) . (lie) 
A 

Thus, to first order in small quantities, the ellipticity is exactly what would be calculated by 
assuming the magnetic field to be frozen at its instantaneous position during its traversal by 
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the laser beam, using the formula for ellipticity calculated from the vacuum birefringence in 
a static magnetic field. 



5. Perturbation method for a rotating, time-dependent magnetic field 

The simple form of the final answer of Eq. (10b) suggests there should be a simpler 
derivation than the wave matching procedure employed in Sec. 3. We show in this section 
that the same answer can be obtained by a perturbation theory approach, which while 
not yielding the structure of the propagation modes in the magnetic field region, has the 
advantage that it extends to the case when the magnetic field strength (that is, the coefficient 
^) is also time dependent. Substituting the general Ansatz 

e^Ei{z,t)B + E2iz,t)z X B , 

(12a) 

b^Bi{z,t)B + B2{z,t)zx B , 

into Eq. (la), and using Eq. (lb) but now allowing ^ to have a t and z dependence, we get 
the coupled equations 

d,Ei{z,t)^-flBi{z,t)-dtB2{z,t) , 

-d,E2{z,t) =-dtBi{z,t) + nB2iz,t) , 

(1 - 60dzBi{z, t) =(1 + 5^)nEi{z, t) + (1 - 2^)dtE2iz, t) - 2dt^E2{z, t) + 6d,^B,{z, t) , 

-(1 - 2^)d,B2{t) =(1 + bmEi{z,t) - (1 - 20nE2{z,t) + 5dtCEi{z,t) - 2d,^B2{z,t) . 

(126) 

We can solve these equations by making a perturbation expansion in powers of ^, by writing 
-E'1,2 = -Ei°2 + + similarly for Bi^2-, with the fields with subscript (0) the incident 
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wave components on rotating bases, 

^(0) ^gia>(z-t) _ ^ 

(13a) 

^(0) ^ _ ^iu,{z-t) _ j^^) ^ 

^(0) ^gi^(.-t) ^Qg^^ _ ^ 

and with the fields with superscript (1) perturbations proportional to ^. The zeroth order 
fields are exact solutions of Eq. (12b) when ^ is set to zero, but fl is kept non-zero. Substi- 
tuting the perturbation expansion into Eq. (12), differentiating dzE^l with respect to z to 
get d'^E[^l, and using Eq. (12b) to eliminate cross derivatives dtdzB[^l in terms of e[^1, we 
get inhomogeneous wave equations for the first order perturbations, 

h ^l^d^Ef^ + l2dtidtEf^ + bdl^Ef^ - 2dzidtEf^ - 2dME? , (136) 

h =4ea,2^f + 2dtidtEf^ - 2dl^Ef^ - QdzidtEf^ - Qd^^Ef^ , 
where we have dropped terms of order (but have kept Q in the phases of the zeroth order 
solution, where it multiplies t, which can be arbitrarily large). Using the Green's function 
G{z, t) for the one-dimensional wave equation with outgoing wave boundary conditions, 

G{z,t)^-\e{t-\z\) , 

2 (13c) 
{dl-dl)G{z,t)=5{z)5{t) , 

we can solve Eq. (13b) to give 

pL poo 

Ei]l{z,t)^ dz' dt'G{z-z',t-t')h,2{z',t') . (14) 

Jo J-oo 

A similar solution can be obtained for the perturbations B^l, giving the general first order 
corrections to the fields, even when the magnetic field strength is t and z dependent. 
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Let us now consider the case when the field strength is only weakly varying in time, 
so that dt^/^ is a small parameter. Then to first order in small quantities, the terms in 
Eq. (13b) involving time derivatives of ^ as well as explicit powers of can be dropped, so 
that /i,2 reduce to 

h =e^^^^-*) cos(^ - nt) [-7^0)^ + 2iujd,i\ , 

(15a) 

/2 =e^^(^-*) sin(^ - VLt) [-A^uj'^ + Uujd.^] . 
When the magnetic field is uniform (although possibly weakly time dependent) in the interval 
< 2; < L, and zero elsewhere, one has d^^ — ^{t)[5{z) — 5{z — L)], with dt^{t) second order 
small. Then substituting Eq. (15a) into Eq. (14), taking 2; = L, so that \z — z'\ = L — z' 
on the support of the integrand, and regarding terms of order Q/cD and flL as small, we can 
evaluate the integrals over z',t' to get the first order fields. For the contribution from the 
terms — (7, 4)^a)^ in the square brackets in Eq. (15a), which have support on the interval 
< z' < L, we use 

/ dz' / dt'9{t -t'-L + ^')e-(-'-*')±^f^t' 

Jo J -00 



r dz''-t 

Jo ^ 



t-L+z' 



~ I _^iu>{L-t)±iQ{t-L+z') 



(156) 



— _ c; , 

while for the contribution from the terms (2,6)^(1)9^^ in the square brackets in Eq. (15a), 
which have support on the field boundaries at 2;' = 0, L, we use 



dz' / dt'e{t -t' -L + z')^{t')[6{z') - 6{z' - L)]e 



iuiiz' -t')±int' 



j diie''^(^-")±^^"[e(M-L)e^^"^-e(M)] (1^^) 
J —00 
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Adding back the zeroth order fields, we get for the fields at z — L, 

Ei\z=L =e*^(^-*) cos(^ - nt){l + liu^L) , 
E2\z=L =e'^("-*) sin(^ - VLt){l + Uoj^L) , 

Zi 

which when transformed back to the fixed bases by 

Eil=L{cosQt,smQt,0) + E2\z=L{-smQt,cosQt,0) = e'^'-''^\X,Y) 

give 

11 3 

X = cos 9{1 + —iOLC) + -^iojL^ cos(2Qf - 9) , 

11 3 
y = sin ^(1 + —iuLC) + -iuL^ sin(2Qt - 9) , 

in agreement with our earlier result of Eq. (10b). 

6. Discussion 

Returning to Eq. (5b), we see that the form of the propagation eigenmodes in the 
magnetic field region is governed by the behavior of the square root term in a±, that is, by 
((20/a;)^ + (3(^/2)^) ^, with distinctly different behaviors in the small f2 regime 2n/uj << 
3^/2, and the large Q regime 2Q/u; » 3^/2. In the former, the refractive indices are 
~ 1 + (7/2)^ and n_ ~ 1 + 2^, as in the case of a non- rotating magnetic field. In the latter, 
the refractive indices are n± ~ l + (ll/4)^±Jl/a;. This change in character of the eigenmodes 
does not show up in the transmitted wave polarization parameters, however, because all 
dependence on the square root cancels between the various terms in Eq. (11a). The simple 
results that we get for the polarization parameters are what one would immediately get by 
assuming that the dependence on ^ and ^2 should be analytic around the origin in these 
parameters, since then the fact that effects of the magnetic field region vanish at ^ = 
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(16a) 



(166) 



(16c) 



implies that there can be no first order term proportional to Q; the rotation dependence 
must first enter at second order through terms proportional to ^fl. Similarly, dependences 
on a small time rate of change of the magnitude of the magnetic field, which are relevant for 
experiments in which a non-rotating but ramping magnetic field is used, must also appear 
only in second order terms in small parameters. 

We note in conclusion that putting in the numbers characterizing the PVLAS exper- 
iment, even though their magnetic field rotates only a fraction of a revolution per second, 
their parameters yield QL = 0.7x 10-^ 2Q/u} ~ 2 x 10-l^ and 3^/2 ~ IQ-^^. So although all 
of these are small, PVLAS is in fact operating in the large fl regime in terms of the behavior 
of the square root and the propagation eigenmodes in the rotating field region. Nonethe- 
less, as we have shown, the formulas they have used for the ellipticity and polarization axis 
orientation of the emerging wave are correct. 

Acknowledgments 

This work was supported in part by the Department of Energy under Grant #DE- 
FG02-90ER40542. I wish to thank Raul Rabadan for bringing the paper of ref [6] to my 
attention, and participants in the "Axions at the Institute for Advanced Study" workshop 
for stimulating discussions. After the initial draft of this paper was posted to the arXiv, 
I benefited from informative email correspondence with Carlo Rizzo and Holger Gies, and 
from helpful referee comments. 

Added note: After this work was completed, we learned of related work by Biswas and 
Melnikov [10], whose results agree with ours. 

18 



References 



[1] Heisenberg W and Euler H (1936) Z. Physik 98 714 
[2] Schwinger J (1951) Phys. Rev. 82 664 
[3] Toll J (1952) unpublished Princeton dissertation 
[4] Adler S L (1971) Ann. Phys. 67 599 

[5] Zavattini E et al. (2006) Phys. Rev. Lett. 96 110406; Bakalov D et al. (1994) Nucl. 
Phys. B Proa. Suppl. 35 180; Bakalov D et al. (1998) Quantum Semiclass. Opt. 10 239; 
Bakalov D et al. (1998) Hyperfine Interact. 114 103. 

[6] Tito Mendonga J, Dias de Deus J and Castelo Ferreira P (2006) Phys. Rev. Lett. 97 
100403 

[7] Tito Mendonga J, Dias de Deus J and Castelo Ferreira P (2006) |arXiv:liep-ph /060931 1 



[8] Problems with the calculation of ref [6] were first pointed out by Holger Gies (unpublished 
memorandum, July, 2006), and problems with ref [6] and ref [7] were noted independently 

19 



by the author (unpubhshed memorandum, October, 2006) 



[9] Born M and Wolf E (1964) Principles of Optics (New YorkiMacmillan); for a concise 

exposition with apphcation to PVLAS, see 

Adler, S L (2006) Pittp:/ /www.sns.ias.edu/^adler/talks.htrn] 



[10] Biswas S and Melnikov K (2006) arXiv: |hep-ph/06ll"515 



20 



